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ABSTRACT 

Context. Stars form in the cold dense cores of interstellar molecular clouds and the detailed knowledge of the spectrum of masses of 
such cores is clearly a key for the understanding of the origin of the IMF. To date, observations have presented somewhat contradictory 
evidence relating to this issue. 

Aims. In this paper we propose to derive the mass function of a complete sample of dense molecular cores in a single cloud employing 
a robust method that uses uses extinction of background starlight to measure core masses and enables the reliable extension of such 
measurements to lower masses than previously possible. 

Methods. We use a map of near-infrared extinction in the nearby Pipe dark cloud to identify the population of dense cores in the cloud 
and measure their masses. 

Results. We identify 159 dense cores and construct the mass function for this population. We present the first robust evidence for a 
departure from a single power-law form in the mass function of a population of cores and find that this mass function is surprisingly 
similar in shape to the stellar IMF but scaled to a higher mass by a factor of about 3. This suggests that the distribution of stellar birth 
masses (IMF) is the direct product of the dense core mass function and a uniform star formation efficiency of 30%±10%, and that the 
stellar IMF may already be fixed during or before the earliest stages of core evolution. These results are consistent with previous dust 
continuum studies which suggested that the IMF directly originates from the core mass function. The typical density of ~10 4 cirr 3 
measured for the dense cores in this cloud suggests that the mass scale that characterizes the dense core mass function may be the 
result of a simple process of thermal (Jeans) fragmentation. 
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1 . Introduction form differs from the form of the stellar IMF in two fundamental 

ways. 

The stellar initial mass function (IMF) is one of the most fun- 

damental distributions in astrophysics. Detailed knowledge of Flrst ' the stellar IMF cannot be described by a single power- 
its functional form and how this form varies in space and time law (scale-free) function. The general consensus of observa- 
is necessary to predict and understand the evolution of stellar tlonal studles 18 mat me stellar IMF 18 best described by a se- 
systems, from clusters to galaxies. Nevertheless, the origin of nes of Power-law segments that overall is qualitatively similar 
the stellar IMF remains one of the major unsolved problems of "Lshape to a log-normal function ( |Miller & Scalo| | 1979t | Kroupa| 
modern astrophysics 120021 ; IMuench et al.1120011) . In particular, the IMF is character- 
Stars form in interstellar molecular clouds and the detailed ^ ed by a broad peak which spans roughly the interval of 0.1-0.6 
knowledge of the spectrum of structure in such clouds is likely M ° ™ d ™\ °^i pidly ° n each Slde ° f * e P f k " Between 06 
a key for the understanding of the origin of the IMF. What and 10 M ° I, , decreasgsma E2wer-law fashion similar to 
is the relation between cloud structure and the IMF? To date, * at originally derived bylSaJpstej dl955|). The fact that the IMF 
observations have presented somewhat contradictory evidence dls P[ a y s a P eak . and . ls . not scale-free is of great significance for 
relating to this issue. Early CO observations suggested that star Ration since it indicates that there is a characteristic mass 
the mass spectrum of molecular cloud st ructure is we ll de- result ing from the star formation process of around 0.2-0.3 M 
scribed by a single pow er-law (|Bljy|l^|Krarner_eIiD[l9 (e.g.,| Larson| | 2005[ ). 

iHeithausen et all 1998b . specifically, dN(m)/dlogm ~ m~ 0J for Second, above about a solar mass the slope of the stellar 

structures (clumps) with a wide range of masses that apparently IMF is steeper (i.e., dN(m)/dlogm ~ m~ 135 ) than the cloud 

extends from a fraction a solar mass to 1000 M Q . This functional mass spectrum. This latter difference implies a significant phys- 

ical difference in the two distributions. For the stellar IMF, the 

Send offprint requests to: J. Alves bulk of the (stellar) mass is tied up in low mass objects while for 

* University of Milan, Department of Physics, via Celoria 16, I- clouds the bulk of the mass is tied up in the most massive objects. 

20133 Milan, Italy A consequence of this difference is that in order to produce the 
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Fig. 1. Dust extinction map of the Pipe nebula molecular complex from iLombardi et aD ([2006). This map was constructed from 
near-infrared observations of about 4 million stars in the background of the complex. Approximately 160 individual cores are 
identified within the cloud and are marked by an open circle proportional to the core radius. Most of these cores appear as distinct, 
well separated entities. 



stellar IMF from the cloud mass spectrum, a transformation must 
take place during the process of star formation. It has been sug- 
gested that outflows generated during the protostellar stages of 
star formation provide a natural feedback t o collapse/infal l lim- 
iting the final mass of a protostellar object dShu et al 1 119871) and 
perhaps providing the mechanism for transforming the form of 
the cl oud mass spectrum into the form of the stellar mass spec- 
trum (lLada & Ladall2003l) . However, it is not clear how such a 
process could produce a peak, or a characteristic scale for stellar 
masses. Moreover, the comparison of the CO core mass function 
and the IMF may not be relevant since CO does not trace the 
dense co mponent of t he molecular gas within which stars actu- 
ally form dLadall 19921) . 

Indeed, a different picture appears to emerge when obser- 
vations of dense gas are considered. Typically only about 10% 
or less of the mass of a star-forming molecular cloud is in the 
form of dense (i.e., n(H2) ~ 10 4 cirT 3 ) gas and this gas ap- 
pears to be organized i nto d iscrete cores within which stars form. 
Tachi hara et all (|2002|) and lOnishi etail (120021) . using C ls O and 
H 13 CO + as tracers of dense gas, suggest that the stellar and 
core mass distributions are similar. Recent observations of dust 
continuum emission originating from such dense cores has en- 
abled the construction of the dense core mass function (DCMF) 
in a number of nearby molecular cloud complexes. For cores 
with masses in excess of -0.5 M the derived mass spectra ap- 



pear to be described by a single power-law, similar to the CO 
cores but with a relatively steep slope (-1.1 t o -1.6, in log mass 
units) similar to the that of the stellar IMF (jMotte et al.Nl 998 ; 



Testi & Sargentlll998t Uohnstone et al.ll2000l l200U iMotte et ; . 
200U iBeuther & Schilkd 120041: IStanke et alll2006l) . Moreover, 



in one example, the p Ophiuchi cloud, the core mass spectrum 
exhibited possible evidence of a flattening or break near a mass 
of ab out 0.5 M , also similar to the stellar IMF dMotte et alj 
1998). In another cloud, NGC 1333, measurements of dust emis- 
sion produced a core mass spectru m between 0.1 and .5 M Q 
with a slope of approximately -0.4 dSandell & Kneel200ll) . sim- 
ilar to that universally derived for less-dense gas traced by CO 
emission in other clouds but also consistent with the apparent 
break and assumed flattening below 0.5 M of the DCMF of the 
Ophiuchi dark cloud mentioned above. The observed similarity 
between the slopes of the DCMF derived from millimeter-wave 
dust emission studies and the slope of the IMF above 0.5 M Q has 
been taken as evidence that the individual dense cores are the 
direct precursors of new stars and moreover that the stellar IMF 
is completely specified by the fragmentation process in molecu- 
lar clouds. In addition, the possible flattening of the DCMF near 
0.5 M implies a high star formation efficiency (SFE) for dense 
gas (about 100%). In this case there is no need for a mass trans- 
formation from the DCMF to the IMF. The characteristic scale 
of stellar mass demanded by the stellar IMF is set by the funda- 
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Fig. 2. Mass function of dense molecular cores plotted as filled 
circles with error b ars. The grey line i s the stellar IMF for the 
Trapezium cluster (iMuench et al .1 12002b . The dashed grey line 
represents the stellar IMF in binned form matching the resolu- 
tion of the data and shifted to higher masses by about a factor of 
4. The dense core mass function is similar in shape to the stel- 
lar IMF function, apart from a uniform star formation efficiency 
factor. 



menta l physics of cloud fragmentation dShu et al" 20041 Larson! 
120051) . Which picture is correct? 

In order to bring new insight to this issue we used an in- 
dependent method of identifying and measuring the masses 
of dense cores. This method uses precise infrared measure- 
ments of dust extinction toward stars background to a molec- 
ular cloud (t h e Nice(r) method. lLada et al.l dl994b : lAlves et all 
d!998l 12001b : ILombardi & Alvesl d2001bT Such measurements 
are free from many of the complications and systematic un- 
certainties that plague molecular-line or dust emission data and 
thus enable robust maps of cloud structure to be constructed. We 
used data from the re cent wide-field extinct ion map of the Pipe 
Nebula constructed by Lombardi et al. (2006), hereafter LAL06, 
using 2MASS data. The Pipe n ebula is a virtually unstudied 
nearby molecular cloud complex (Onishi et al. ( 1999), LAL06), 
at a distance of about 130 pc and with a total mass of ~ 10 4 M Q . 

2. Observations 

The details of the extinction study are described in LAL06. 
Briefly, this molecular complex was selected because 1) this is 
one of the closest to Earth complex of this size and mass, 2) 
it is particularly well positioned along a relatively clean line of 
sight to the rich star field of the Galactic bulge, which given 
the close distance of the Pipe nebula allowed us to achieve spa- 
tial resolutions of ~ 0.03 pc, or about 3 times smaller the typ- 
ical dense core size, and 3) it exhibits very low levels of star 
formation suggesting that its dense cores likely represent a fair 
sample of the initial conditions of star formation. LAL06 ap- 
plied a 3-band (1.25 //m, 1.65 fim, 2.2 /im) optimized version 



of this method, the Nicer method (ILombardi & A lves 2003), to 
about 4 million stars background to the Pipe nebula complex to 
construct a 6° x 8° dust column density map of this complex, 
presented in FigureQ] Because of the high dynamic range in col- 
umn density achieved by this map (3<x ~ 0.5 < Ay < 24 mag or 
10 21 < N < 5 x 10 22 cm -2 ), cores are easily visually identified 
as high contrast peaks embedded on rather smooth but variable 
background (see Figure [TJ. Unfortunately, because of the high 
dynamic range and variable background, traditional source ex- 
traction algorithms based on thresholding fail to identify these 
objects in a coherent way. An alternative approach is to extract 
cores based on their sizes, using a multi-scale algorithm. For this 
study we used a algorithm developed by Vandame (2006, private 
comm.), which, in brief, uses the wavelet transform of the image 
to first identify and then reconstruct the dense corefl This step 
defines the projected core boundaries. 

Masses are derived by integrating the extinction map over the 
area of each core and multiplying by the standard gas-to-dust ra- 
tio. The final Pipe core sample has 159 objects with effective di- 
ameters between 0. 1 and 0.4 pc (median size is 0. 18 pc) and peak 
extinctions that range from 3.0 to 24.3 visual magnitudes (mean 
extinction is 8.4 magnitudes). The derived core masses range be- 
tween 0.5 to 28 M . The assessment of sample completeness is 
non trivial because of the variable background. Nevertheless, the 
completeness should not be dominated by confusion but sensi- 
tivity, as the mean separation between cores, even in the clus- 
tered regions, is well above the resolution of the map. We esti- 
mate, conservatively, that the sample is 90% complete at about 
1 M . The mean diameter of a 1 M Q core is ~0.2 pc, i.e., about 
six times the resolution of the map. 

3. Results: The Dense Core Mass Function (DCMF) 

The dense core mass function we derive from the above obser- 
vations is presented in Figure [20. For comparis on we plot the 
Trape zium cluster stellar IMF as a grey solid line (M uench et al.l 
2002). This IMF consists of 3 power law segments with breaks 
and 0.6M G and O.12M We find that the DCMF for the Pipe 
Nebula is surprisingly similar in shape to the stellar IMF. 
Specifically, bot h distribution s are characterized by a Salpeter- 
like power-law dSalpeterll 1955b that rises with decreasing mass 
until reaching a distinct break point, this is then followed by a 
broad peak. Although similar in shape, the stellar and core mass 
functions are characterized by decidedly different mass scales. 
The grey dashed line in Figure [2] is not a fit to the data but sim- 
ply the stellar IMF in binned form matched to the resolution of 
the data, and shifted by a factor of 4 to the higher masses. The 



1 Object identification in wavelet space: for a given scale i, structures 
are isolated with classical thresholding at 3<r, with cr, being the noise 
amplitude at scale i. A structure at scale i is connected with a structure 
at scale i + 1 if its local maxima drops in the structure at scale i + 1. 
The size scales considered were 2', 4', and 8' (0.08 pc, 0.15 pc, and 
0.30 pc). One then builds a 3D distribution of significant structures (x, 
y, and i). The algorithm developed by Vandame (2006, private comm.) 
offers rules that split the 3D distribution into independent trees corre- 
sponding to one core and its corresponding hierarchical details. Object 
reconstruction: This same algorithm performs a complex iterative re- 
construction of the cores following the trees defined in the previous 
step. The final "cores only" image is validated by subtracting it from 
the extinction map, effectively creating a smooth image of the variable 
background. 

2 The full dataset is available in electronic form at the CDS 
via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat ?J/A+ A/ 
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break from the Salpeter-like slope seems to occur between 2 and 
3 M H for the cloud core s instead of the 0.6 M for the stellar 
case dMuench et al.ll2002l) . 

The DCMF in Figure [2] likely suffers from two sources of 
uncertainty: 1) the individual core masses are likely upper lim- 
its to the true values since we made no corrections for the local 
extended background on which most cores are embedded, and 
2) the particular binning we chose may not produce the most 
accurate representation of the underlying mass distribution. To 
address point 1) we estimated a lower limit to the true core 
masses (and the DCMF) by subtracting from each core its lo- 
cal background. We then constructed a background subtracted 
DCMF and found that the form of the mass distribution was 
maintained although the new distribution was shifted towards 
the low masses. To address point 2) we compute a better func- 
tional representation for the DCMF using a probability density 
estimator that does not require data binning. We used a gaus- 
sian k ernel estimator w ith a window width of 0.15 in units of log 
mass (Hlverman 1981). 

We present in Figure [3] the probability density function for 
the core masses (background subtracted). As already suggested 
in Figure 2, the core mass distribution seems to be characterized 
by two power-laws and a well defined break point at ~ 2 M . 
Fo r comparison w e also present th e field sta r IMF determined 
by iKroupal j200l) (solid grey) and IChabrierf J2003h IMF (dot- 
ted gr ey), and the IMF for the young Trapezium cluster (dashed 
grey) dMuench et al.ll2.Q02D already present in Figureff] All these 
stellar IMFs were shifted to higher masses by an average scale 
factor of about 3. The point of this comparison is to simply il- 
lustrate the overall similarity between different stellar IMFs from 
different environments and constructed in different ways, and the 
DCMF of the Pipe Nebula. Considering the overall uncertain- 
ties, apparent differences between the different distributions are 
likely not significant. While Figure [3] presents the background 
subtracted core sample, the non background subtracted sample 
is virtually identical in shape, although shifted towards higher 



4. Discussion and Conclusions 

The close similarity in shape between the DCMF and the IMF 
supports a general concept of a 1-to-l correspondence between 
the individual dense cores and soon to be formed stars. In this re- 
spect our observations are consistent with and appear to confirm 
the results of dust continuum surveys in other clouds. However, 
our study provides the first robust evidence for a departure from 
a single power law in the core mass function. We find that the 
location of the break indicates that there is a factor of about 3 
difference in mass scale between the two distributions, which 
implies that a uniform SFE, that we estimate to be 30%+ 10%, 
will likely characterize the star formation in these dense cores, 
across the entire span of (stellar) mass. This efficiency is very 
similar to the maxim um estimated for young embedded clusters 
dLada & Ladall2003l) . 

It is well known that the generation of a powerful outflow is 
a nat ural conse quence of the formation of a protostar in a dense 
core (Lada 1985). It has been long suspected that such outflows 
could be very disruptive to surrounding material and could play 
an importan t role in limiting the mass which can accre t e onto 
a protostar dShu et alJll987lll99ll) . iMatzner & McKed d2000l) 
have theoretically investigated the disruption of dense cores by 
protostellar outflows and have concluded that outflow disruption 
of individual cores can produce SFEs in the range of 25-75%. 
Moreover they found that SFE is independent of stellar mass, 
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Fig. 3. Probability density function of core masses (black) for the 
background subtracted sample. The grey region indicates sam- 
ple incompleteness while the two thin dashed line indicates the 
90% confidence limits. The core mass distribution seems to be 
characterized by two power-laws and a well defined bre ak point 
at aro und 2 M . Also pl otted are t he fie ld star IMFs of Kroupa 
d200ll) (solid grey) a nd IChabrier! d2003l) (dotted grey), and the 
Mue nch et alj d2002l) (dashed grey) IMF all scaled up by a fac- 
tor of ~ 3 in mass. 



therefore the stellar and core mass fu nctions should b e of simi- 
lar overall functional form. Recently, IShu et all d2004l) have de- 
veloped a self-consistent analytical theory for the origin of the 
stellar IMF incorporating the feedback from outflows. In their 
picture, cores are formed from magnetized turbulent gas and are 
disrupted by outflows with a resulting SFE predicted to be « 
30%. The close correspondence of this prediction with our esti- 
mated SFE for the Pipe cores may implicate outflows as the key 
mechanism for core disrupti on and for sett ing the final masses 
of the protostellar objects. IShu et al.l d2004) also predicted both 
a slope for the power-law portion of the core mass function and 
a location of the break from the power-law form that are very 
similar to those we measure in the Pipe Nebula. 

A number of other theoretical calculations have also 
produced DCMFs s imilar in shape to that f ound here 
dPadoan & Nordlundl d2002l) ; iBonnell etail d2006l) : lElmegreenl 
(2006)). In particular, numerical simulations of turbulent clouds 
predict DCMF shapes similar to the stellar IMF. However it is 
difficult to compare these calculations in more detail to the de- 
rived Pipe DCMF because, unlike the IShu et all d2004h model, 
the mass scales for these DCMFs are presently unconstrained 
and appear to be dependent on model parameters, such as the 
numerical resolution of the simulations as well as the Mach num- 
ber o f turbulence wh ic h are difficult to quantify observationall y 
(e.g., iKlessenl d2001h ; iBallesteros-Paredes & Mac Low! d2002l) : 
iGammie et al.l d2003l) ). Nonetheless we expect that with appro- 
priate modification of model parameters such models could also 
be made consistent with the observations. 
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These considerations suggest that apart from a uniform effi- 
ciency factor, which is likely fixed by the generation of an out- 
flow, the birth mass of a star could to be completely predeter- 
mined by the mass of the dense core in which it is born. 

The existence of a characteristic mass in the stellar IMF does 
suggests a characteristic mass scale is produced by the physi- 
cal process of star formation. The physical origin of this mass 
scale has not been clear although it has been long known that 
thermal (J eans) fragmentation can produce a mass scale for star 
formation (lLarsonll 19851 1 1986b . For a pressure bounded core the 
appropriate mass scale is the critical Bonnor-Ebert mass, 

M BE * 1.15 X («/10V /2 x (7710) 3/2 (1) 

where T is the gas temperature and n its volume density. The 
characteristic mass scale suggested by the stellar IMF is about 
0.3 M . For a typical dark cloud of temperature of 10 K, the cor- 
responding density scale would need to be about 1.5 x 10 5 cm" 3 . 
This is significantly higher than the measured average mean den- 
sity of the Pipe cores of 8xl0 3 cm 4 . However, if the character- 
istic stellar mass scale is set by the characteristic mass of the 
DCMF, which is a factor of ~ 3.3 higher than that of the stellar 
IMF, then the corresponding critical density is ~ 1.3 x 10 4 cirT 3 
much closer to that observed for the Pipe dense cores. This cor- 
respondence suggests that the fragmentation scale for the DCMF 
and ultimately the stellar IMF may be set by a process of ther- 
mal (Jeans-like) fragmentation, perhaps modified somewhat by 
ambipolar diffusion and/or turbulence. Interestingly, from the re- 
cent analysis of Spitzer observatio ns of the spatial distr ibution 
of protostellar objects in NGC 2264, lTeixeira et all ((2006) found 
evidence for Jeans like fragmentation to be occurring in that star 
forming region. 

We note that the mean densities of the cores in our study 
range between 5xl0 3 and 2xl0 4 cirT 3 (for the background sub- 
tracted case). This range is lower than that (3xl0 4 ~ 5 cm -3 ) 
usually associated with typical star forming cores and even 
the st arless cores typically traced by dust continuum observa - 
tions dMotte et alJl998HJohnstone et al.ll200UlKirk et al.ll2005l) . 
Because extinction observations are more sensitive to lower col- 
umn densities we are able to probe lower mass and density 
regimes than those probed by current dust continuum studies. 
If the Pipe Nebula cores are to form stars, the relatively low 
mean densities they currently possess may indicate that they are 
in a very early stage of development. This notion is somewhat 
strengthened by the absence of significant star formation activity 
in the Pipe nebula complex. Apparently, the functional form of 
the stellar IMF may be fixed during or before the earliest stages 
of core evolution. 

The results of this letter suggest that the distribution of stellar 
birth masses over essentially the entire range of stellar mass is 
the product of the dense core mass function and a uniform SFE 
of 30%+10%. Hence, the question of the origin of the IMF is 
the question of the origin of the DCMF, itself a result of the 
physical process of cloud fragmentation/coalescence. The mass 
scale characterizing the DCMF may result from a simple process 
of thermal (Jeans-like) fragmentation in the molecular gas. 
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